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Abstract 

Thin-shell AdS-Vaidya spacetimes can be considered as holographic mod- 
els of the thermalization process in strongly-coupled conformal field theories 
following a rapid injection of energy from an external source. While the 
expected thermalization time is the inverse temperature, Bhattacharyya and 
Minwalla have pointed out that bulk causality implies that expectation values 
of local field-theory observables actually take on their thermal values imme- 
diately following the injection. In this paper we study two-point functions, 
for which the causality argument does not apply. Specifically, we study the 
Brownian motion of a "quark" represented by a string stretching from the 
boundary to the horizon of an AdS3-Vaidya spacetime. Surprisingly, we find 
that the two-point function also thermalizes instantly. Since Brownian mo- 
tion is a ji effect, our result shows that, at least in certain cases, the rapid 
thermalization property of holographic plasmas persists beyond leading order 
in i 
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1 Introduction 

One of the fundamental problems in physics is to understand the dynamics of a ther- 
mal system out of equilibrium. The time-dependence nature of the problem makes 
it very hard to explore. An interesting dynamical system which has prompted a lot 
of interest is the quark-gluon plasma (QGP) produced at RHIC experiment [1,2]. 
Modeling based on near-ideal hydrodynamics suggests very rapid local thermaliza- 
tion after the collision [3] . The process of local equilibration from an initial far from 
equilibrium state needs to be understood. Since QGP is a strongly coupled non- 
Abelian plasma, AdS/CFT correspondence [4-6] provides a powerful tool to explore 
the dynamics of this system. 

In thermal field theories stable stationary solutions that describe the systems in 
equilibrium reveal the late-time behaviour of the system. In gauge/gravity duality 
this corresponds to a stationary black hole in the dual geometry. An interesting 
question is to explain how the system goes back to equilibrium after being slightly 
perturbed. From field theory point of view this question has been investigated in two 
regimes that have some overlap in the long wave-length limit. One is linear response 
regime which explores the system with very small fluctuations compared to the mi- 
croscopic scale of the system. In this regime the evolution of the system is explained 
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by linear response theory and fluctuation-dissipation theorem [7] . This calculation is 
based on computing the correlation functions in the field theory. AdS/CFT provides 
us with an efficient tool to obtain the correlation function of the operators in the 
CFT using the bulk dual. The calculation of the retarded correlation functions in 
AdS / CFT was first initiated in [8] . More references on this can be found in the very 
nice review [9]. From the statistical point of view the fluctuation-dissipation theo- 
rem has been proved in the AdS/CFT context in the papers by de Boer, Hubeny, 
Rangamani, Shigemori [10] and Son, Teaney [11]. They have studied the Brownian 
motion of the quark on the boundary of BTZ black hole using the bulk information. 

The linear response of the system in field theory translates into the behaviour 
of the linearized fluctuations of the bulk fields on AdS black hole background. The 
first step in this direction is to study the connection between quasi-normal modes of 
the black hole and the rate at which the fluctuations return back to equilibrium [12]. 

The other regime which has a more coarse-grained view of the problem is the 
hydrodynamic or more generally fluid dynamic regime. This describes the long 
wave length or low energy physics in which the local fluid variables vary very slowly 
compared to the microscopic scale in a way that preserves the local thermal equi- 
librium of the system every where. In the gauge/gravity picture this regime has 
been explored in the interesting works on fluid/gravity duality which was initiated 
by [13]. 

Now one can ask what happens when the system is far from equilibrium which 
does not sit in any of the regimes above. It has been argued in the literature that a 
system out of equilibrium corresponds to black hole formation in the bulk. Recently 
this problem has been addressed in interesting work by Chesler and Yaffe [14] and 
also the work by Bhattacharyya and Minwalla [15]. The former has numerically 
studied the thermalization in M = 4 SYM by triggering the gravitational wave 
and perturbing the boundary stress tensor. The latter has analytically explained 
the formation of the black hole in the bulk using both minimally coupled scalar 
wave and graviton perturbation which corresponds to the injection of energy on the 
boundary. They generally argue, based on the bulk causality, that in the process 
of the black hole or black brane formation the local operators on the boundary 
thermalize instantly while the non-local operators such as Wilson loops thermalize 
in longer time which depends on the scale of the loop. We will base our set-up on 
the result that a rapid injection of energy on the boundary is dual to a thin-shell 
Vaidya spacetime [15]. This seems a reasonable toy model to study thermalization 
on the boundary. 

Motivated by the work [15], in this paper we are taking the first step towards ad- 
dressing thermalization of non-local operators in the holographic field theory where 
the bulk metric is planar Vaidya spacetime corresponding to the planar black hole 
formation. An interesting set-up to address this question is to study the Brownian 
motion of the quark on the boundary, similar to [10]. The quark is the end point 
of a string stretched between the boundary and the horizon. Its fluctuations are 
described by the Hawking radiation of the black hole [10, 11] (The basic idea that a 



2 



string attached to a black hole fluctuates due to the Hawking radiation is discussed 
in [16, 17]). As a candidate for non-local operators we will investigate the two point 
function of the scalar field explaining the fluctuations of the quark. Due to the 
causal structure of the Vaidya spacetime, retarded Green's function will thermalize 
instantly. By this we mean that it equals its thermal value when both arguments 
are taken to be at any time (no matter how soon) after the injection of energy. 
But one might initially think that the Hadamard function of the field which is the 
expectation value of the anti-commutator will behave differently. This is due to the 
fact that in order to calculate the Hadamard function one needs to know the state 
of the thermal vacuum which encodes the information about the Hawking radiation 
of the black hole. In general one expects that it takes some time for the Hawking 
radiation to turn on. 

One of the obstacles that we face is how to calculate the Hawking radiation at 
finite time where the bulk spacetime is the formation of a black hole. The usual 
method used in the literature is to compute the Bogoliubov coefficients. In general 
this is a very hard calculation and has been done only in high frequency or late time 
limit where the spacetime is stationary. This calculation relies on using geometric 
optics approximation. But in this paper we would like to address the on-set of 
Hawking radiation which means calculating a time dependent Hawking radiation 
that in the late time limit reduces to the constant result for stationary black holes. 
In this regard an elegant work has been done by Callan, Giddings, Harvey and 
Strominger (CGHS) [18] where they use the trace anomaly [19] to calculate the time 
dependent Hawking radiation in the formation of 2-dimensional linear dilaton black 
hole. This method can not be used in our system due to the lack of diffeomorphism- 
invariance in the 2-dimensional worldsheet calculation. 

In this paper we introduce a new way to overcome this problem and obtain the 
Hadamard function of the scalar field in a general time-dependent background. This 
is a fairly simple method but one that has not been used in the literature before, as 
far as we are aware. The idea is to calculate the Hadamard function of the fields in 
the vacuum or the spacetime before the formation of the black hole. And then, using 
the fact that the Hadamard function solves the equation of motion in both variables, 
we propagate it from the shock wave to the black hole region using the black hole 
retarded Green's function. We thus replace the problem of solving the wave equation 
an infinite number of times, in order to find a complete set of solutions and compute 
the Bogoliubov coefficients, with the problem of solving the wave equation just once 
but in two variables, which is much more tractable in a general time- dependent 
background (especially a simple one such as thin-shell Vaidya). We will introduce 
and check this method with the known result of CGHS in subsection 4.1. 

The set-up of the paper is the planar BTZ black hole formation in the bulk. We 
will show, using the method mentioned in the previous paragraph, that the two- 
point function or more specifically the Hadamard function of the fluctuation modes 
which is a non-local operator thermalizes instantly after the shock wave (subsection 
4.2). This is an interesting result which seems to be at odds with the arguments 
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of [15]. A simple observation that makes this result less surprising is the general 
discussion of a Rindler observer. If an observer who is at rest in flat space and 
starts to uniformly accelerate, calculates the two-point function of the fields on its 
worldline he will measure a thermal two-point function. 

We will also see that using the method of propagation, the same result of instan- 
taneous thermalization of the Hadamard function is obtained in the 2-dimensional 
set-up of AdS2 black hole formation. This is not a new result and has been previ- 
ously observed in an interesting paper by Spradlin and Strominger [20]. This result 
applies to the transverse fluctuations of the string in unwarped directions, such as 
the S 3 or M 4 directions in AdS 3 x S 3 x M 4 where M 4 can be T 4 or K3. 

To clarify more the result of this paper we mention that if one wants to use 
holographic models to understand (even qualitatively) the behaviour of the real 
quark-gluon plasma, such as its apparent rapid thermalization, one has to under- 
stand the effect, which are potentially very important because in QCD N is just 
3. Whereas previous works studying thermalization in holographic plasmas have 
considered the classical dynamics of the bulk spacetime, here we are considering a 
quantum effect, which from the boundary point of view is a 4 effect. While the 
particular model we study is not a realistic model of the QGP, it is suggestive that, 
in our toy model, the rapid thermalization survives at the level of effect. 

An immediate question that comes to mind is to study the same set-up but in 
higher dimensional Vaidya spacetimes to see if one can observe the same result of 
instantaneous thermalization. Or this result is particular to 3-dimensional setting. It 
can be asked why after doing a complicated calculation for the Hadamard function in 
AdS3 Vaidya background we get such a simple answer. Is there any way to organize 
the calculations that makes the behaviour of the Hadamard function clear. One can 
also look at the Hadamard function of the scalar field on the Vaidya background 
instead of the 2-dimensional worldsheet calculation. This can be done in the 3- 
dimensional set-up of this paper or other higher dimensional black holes. One may 
also ask what happens to the bulk Hadamard function after the shock wave. In 
the AdS2 black hole formation the bulk Hadamard function thermalizes instantly 
too. But we have not done this calculation in the BTZ formation set-up. Another 
interesting approach to thermalization is to study other non-local operators such as 
Wilson loops. From field theory point of view one might also ask if this observation 
only applies to holographic strongly coupled plasmas or in any strongly coupled field 
theory this rapid thermalization happens. All of these paths are very interesting to 
investigate. 

2 The set-up of the problem 

In this section we will explain the set-up of the problem. We would like to study 
the Brownian motion of a quark on the boundary where the bulk dual is the 3- 
dimensional planar Vaidya background. Regarding the result of [15] this can be 
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a good toy model to study thermalization of non-local operators. The quark on 
the boundary is in fact the end point of a string which is stretched between the 
boundary and the horizon. Therefore the fluctuations of the quark are described by 
the fluctuations of the transverse modes of the string. These fluctuations originate 
from the Hawking radiation of the black hole which propagates along the string 
[10,11]. In fact it has been shown that the modes on the string worldsheet which 
attaches to the horizon are thermally excited by the black body radiation of the 
black hole [16, 17]. In this section we will closely follow the notation of [10]. 
The background metric can be written as 

ds 2 = g^{x)dx^dx v + G(x)dX 2 , (2.1) 

where x M = t, r. We would like to work in the static gauge in which the worldsheet 
coordinates are identified with the spacetime coordinates. Therefore one can assume 
that the string is streched along the r direction and its small fluctuations are in the 
transverse direction X. We consider g s , the string coupling, to be very small and 
work in the limit where the motion of string is described in the probe approxima- 
tion which means the backreaction of the string on the background geometry can 
be ignored. Thus the string cr-model action in the absence of the B-field can be 
expanded to quadratic order in X as 

Sng = J d 2 x V-d^tj^u , (2.2) 

« \ Sx V=99^Gd^Xd v X , 

where 7 M „ is the induced metric on the worldsheet. With this set-up the Brownian 
motion of a quark on the boundary can be addressed by doing a 2-dimensional 
worldsheet calculation with a Vaidya metric. 

The bulk spacetime is described by a 3-dim Vaidya metric which can be written 

as 

2 2 

ds 2 = -C- - MQ{v))dv 2 + 2dvdr + T -dX 2 . (2.3) 
V I 2 

where M(v) = MQ(v) is a step function 

M(v) = 0, v < , 

M(v) — M, v > , (2.4) 

and M is the mass of the black hole after formation. The metric has been written 
in ingoing Eddington-Finkelstein coordinate, (v,r). The advantage of writing the 
metric in this coordinate system is that it is non-singular on the horizon and the 
coordinates (w,r)are well-defined all over the space, v = t + r* is the ingoing null 
coordinate which ranges over (—00, 00) and < r < 00. r* is called the tortoise 
coordinate and is defined as 

dr* Fg^ 



dr V 9tt 
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Figure 1: Penrose diagram of the BTZ black hole formation. The dashed line is the event 
horizon and the arrow shows the shock wave which denotes the collapse of null matter. 



Therefore we have 

AdS 3 



> 



r 

BTZ : r* = — ln(^^) , (2.6) 

where r# = \/Ml is the black hole horizon radius and — oo < r* < 0. 

The metric (2.3) explains the formation of a black hole with mass M by the 
collapse of null matter (shock wave) at v — 0. The shock wave is stretched in the X 
direction. One can see due to the r 2 term in the third component of the metric, the 
shell shrinks as r — > 0. The bulk spacetime is AdS3 before v = and non-rotating 
planar BTZ black hole after it. By planar BTZ black hole we mean IgI, instead 
of normal definition for BTZ which is X = lip with if = if + 2tt. 
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Basically we would like to work with the null coordinates (v, u) where u = 
t — r* is the outgoing null coordinate and ranges over — oo to oo. But it is not 
continuous across the shock wave and we have to find the relation between outgoing 
null coordinate in AdS 3 and BTZ. Let's define u as a function of u{v,r). Regarding 
the fact that we want to write the metric in null coordinates, we have ds 2 = when 
du = 0. Therefore in order to calculate r as a function of v and u we have to solve 
the differential equation 

(IT 1 

2---(r 2 -^9(,)) = 0. (2.7) 

This equation can be solved for the two regions before and after the shock wave and 
we get 

-2/ 2 

v < , 



v + h(u) 



1 + e ~P 

r = -r H ,. ; v > . 2.8 

I _ e —fi-(v+f a (u)) 

Since r is the well-defined coordinate in all the space-time so it should be continuous 
across the shock wave at v = 0. Imposing this condition on the above result we get 

We have the freedom to choose fi(u) and ^(w) in terms of outgoing null coordinate 
in AdS3 and BTZ spacetimes. If we choose fi(u) = —uj and f^u) = —u where u 
and u are the outgoing null coordinate in AdS 3 and BTZ, respectively, we get 



2/ 2 r 

u = — tanh(-§«) . (2.10) 
r H 2/ 2 



Therefore at the horizon of the black hole which corresponds to u = oo we have 



2/ 2 

uj = —. 

rH 



The Penrose diagram of the spacetime is shown in figure (1). The dashed line is 
event horizon and dash-dotted line is the apparent horizon. The role of the apparent 
horizon in non-equilibrium systems has been mentioned in [14]. 

In the following subsections we explain the mode expansions in each space, sep- 
arately. 



2.1 String Mode Expansion in AdS3 

Since we would like to study the behaviour of the correlation function of the fluctu- 
ations of the quark on the boundary we need to obtain the mode expansion of the 
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transverse mode of the string in the two regions of the Vaidya spacetime. We start 
with AdS3 metric in the Poincare patch which is 

ds 2 = -^dt 2 + -dr 2 + ^dX 2 . (2.11) 

The string is stretched along r direction and the worldsheet coordinates are assumed 
to be (t,r) (static gauge). The 2-dimensional worldsheet equation of motion is 

' dJv=9Gg^d„x)=0, (2.12) 



y~9 

where compared to (2.1) we have 

,2 

G(x)dX 2 = T -dX 2 . (2.13) 
l z 

Note that in a higher dimensional spacetime the equation of motion could be written 

as 

1 -A [V=9 Gu ST d v X J \ = . (2.14) 



V~9 

Using (2.11) the equation of motion reduces to 

1 

r 

This equation can be even more simplified if we assume 



d 2 X + r* 2 d T *{^d r *X) = . (2.15) 



X(t,r) = ^±. (2.16) 

Therefore we get 

-d 2 <P + dU-^<P = . (2.17) 

One can see that this equation is similar to the Schrodinger equation with the 
potential 

_9 —It 1 

= - 5 = - ir . (2.i8) 

It is infinite at r — > oo which is due to the fact that the space-time is asymptotically 
AdS 3 . 

Let's assume that the function X(t, r) is separable. Due to the fact that the 
metric is time-independent we can write the classical solution to the wave equation 
for the string modes in the form 

« w (t,r) = e- iwt / w (r) . (2.19) 
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where 00 > 0. Therefore the string modes are given by the function 



u u (t,r) = A 



(r + iul 2 



+ B 



[r — iul 2 



—iuit 



(2.20) 



where A and B are constants and should be fixed by the boundary conditions. In 
order to have a finite mass for the external particle on the boundary in AdS/CFT 
context, we need to impose a UV cut-off near the boundary. We choose it as r c >> 
1. In order to fix the constant B in the wave equation we impose the Neumann 
boundary condition at the cut-off 



drU^r, i)|, 







Therefore we get 



B = -e 



2iu>r* 



(2.21) 
(2.22) 



which is a pure phase. 

The next step is to normalize the modes in order to fix the free constant A. This 
can be done by imposing the normalization conditions 



(u w ,tv) = -«,<,,) = 5(u - J) ; (u w ,iO) = = 

The definition of the inner product is given by 



(f,9h 



2na' 



dxVhn»G IJ (f I d li g J ' - fy/y*; 



(2.23) 



(2.24) 



where S is a Cauchy surface, is the future pointing unit normal vector to E and 
hij is the induced metric on E. The Cauchy surface is given by S = t — t which is 
a constant t hypersurface. Therefore the definitions for the normalization condition 
reduces to 



i f t 

{u^Uu') = -7^7 J dr -^grriu^dtul, - u^dtuj) , (2.25) 

(2.26) 



which leads to 

A- 1 

col V 2w 

For the future reference we mention that in case we consider the modes to be slightly 
massive (//—>■ 0) we get 

1 [of 

A=^-J—, (2.27) 
where u = \/k 2 + fi 2 . Therefore the mode expansion is given by 
1 



u k (t,r) = 



kl 



2a' 



to 



Skr* 



sinA;(r* — r*) — kr* cosk(r* — r*) 



-ioot 



(2.28) 
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In the limit where we set r* — > which means we move the boundary to infinity the 
solution reduces to 



u k (t,r) 



1 2a' 



sin/cr* — kr* coskr* 



-iwt 



(2.29) 



Now that we have solved the classical wave equation for the transverse fluctu- 
ations of the string in AdS3 background, we can write the quantum field for the 
scalar modes as 

/■oo 

X(r,t)= dk[a k u k (t,r) + a{ul(t,r)] , (2.30) 
Jo 

in terms of annihilation and creation operators which satisfy 

[a k , ay] = [4, a\] = ; [a k , a{,] = 8(k - k') . 



(2.31) 



2.2 String Mode Expansion in BTZ 

The same calculation can be done for the BTZ background 



ds' 



2 ty. 2 72 ,y. 2 



r 2 



H 



The equation of motion for the string fluctuations gets the form 



(2.32) 



-d 2 t X+ — d r *(r 2 d r *X) = 



which gets more simplified to 



I 4 sinh 2 (^) 



. 



(2.33) 



(2.34) 



where is defined in (2.16). Note that similar to the AdS3 case we can read the 
potential which is given by the third term in the above equation. It vanishes at the 
horizon and diverges at infinity, as we expect. 

The solution to the equation of motion is the sum over the ingoing and outgoing 
modes [10] 



A 



r H r + iool 2 ( r — -• 



r H + iujl 2 



r + r H 



" B- TH 



ru — iul 2 



r — iul 2 ( r — Tn 



r + r H 



-iwt 



(2.35) 

Similar to what we did in the AdS3 case, we impose the Neumann boundary condi- 
tion at the UV cut-off and we get 



B 



r H - iujl 2 r 2 H + iul 2 r c 2 - , 

6 

rn + iool 2 r 2 H — iul 2 r c 
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(2.36) 



We also use the normalization condition to fix the coefficient A in the mode expan- 
sion. After some calculations we get 



.4 



/ 



a' 



r H V 2a; 



(2.37) 



In the future calculations we might need to regularize the theory by considering 
the modes to be slightly massive (/i — > 0). Therefore the mode expansion gets the 
form 



Uk 



2a' 



„ikr* 



u r(rH + ikl 2 )(r 2 H — ikl 2 rc) 
r 2 H r + k 2 l r c ) cos k(r* — r*) — kl 2 {r 2 H — rcr) sin k(r* — r* 



(2.38) 



-iujt 



In the limit where we set r* — > which means we move the boundary to infinity the 
solution reduces to 



u k = I 



2a' 



1 



u r{rn + ikl 2 ) 



(kl coskr* + r smkr*)e 



*\„—iuit 



(2.39) 



Similar to the AdS3 case the quantum field for the scalar modes we obtained 
classically is 

/•oo 

X{r, t)= dk [b k u k {t, r) + b\ul (t, r)] . (2.40) 
Jo 

Note that the creation and annihilation operators are not the same as the ones in 
AdS3 space due to the formation of the black hole and the different vacuum state. 
But they satisfy the same relations as (2.31). 



3 Green's Function 

To explain the Brownian motion of the quark on the boundary we need to calculate 
the two-point function which describes the correlation between the two points. Our 
background is a Vaidya spacetime and due to its time dependence, this calculation 
can not be done explicitly. In the following section we will introduce a way to 
do this. But before that we need to obtain the Green's function and Hadamard 
function of the two spaces, AdS3 and BTZ, individually. The commutator (and the 
Green's function which is obtained from the commutator of the fluctuation modes) 
does not depend on the state, because in a free theory it is a c-number. This can be 
easily seen from equation (3.1). But the Hadamard function which is defined as the 
expectation value of the anti-commutator of the fluctuation modes, (3.11), depends 
on the state of the vacuum whether it is thermal. In the AdS3 case this is not a 
problem since the vacuum is not thermal. But in the BTZ background we should 
use the fact that the vacuum is thermal and as we will see the calculation relies on 
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the black body radiation of the black hole. The Hadamard function of the modes 
in BTZ background has been obtained in the limit r* — > 0. These calculations will 
be used later to get the Hadamard function of the modes on the boundary in the 
Vaidya background. 



3.1 AdS 3 Two-Point Functions 

In the previous section we studied the mode expansion in AdS 3 space, equation 
(2.30). We can use it to obtain the Green's function (commutator of the modes 
which is also called Pauli- Jordan or Schwinger function) and the Hadamard function 
of the modes. 



3.1.1 Green's Function 

In general the commutator of the modes is not state-dependent and for the mode 
expansion (2.30) is given by 

POD 

(0\[X(t,r),X(t',r')]\0) = 2i / dk lm(u k (t, r)u* k (t', r')) . (3.1) 

Jo 

Using the solution to the wave equation given in (2.28) one can see that the integrals 
are all finite and we don't need to regularize them. So we consider fi = and 
therefore k = u in (2.28). We can divide the causal structure of our space into 
different regions which are shown in figure (2). The regions are explained by 



I : 


# 

r - 


f r*' - 2r* > At , 




II : 


i * 
\r 


- r*'\ > At > r* + r*' - 


2r c * , 


III : 


* 

r • 


-r*' > \ At\ ; r* + r*' - 


-2r* c 


IV : 


i * 
r 


-r*'\ > -At > r* +r*' 


-2r* c 


V : 


* 

r - 


f r*' - 2r* > -At . 





(3.2) 



Therefore the commutators in different regions of the figure (2) are given by 
/ : ( |[X(t,r),X(t',r')]|0) = ^r* 2 



// : ( |[X(t,r),X(t',r')]|0) = ^(r* 2 + r* /2 -At(4r: + At)) , 

III : (0\[X(t,r),X(t',r')}\0) = — ^r* c At , 

IV : (0\[X(t,r),X(t',r')}\0) = ^^(r* 2 + r*' 2 - At(-4r c * + At)) 
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Figure 2: Different regions of the light cone. The dashed line is r* = r* 



V 



(0\[X(t,r),X(t'y)]\0) 



2ia'n 



(3.3) 



r 

c 



where we have defined At = t — t'. Causality forces us to have zero commutator 
(Green's function) in region III, outside the light cone. But the above result is not 
zero in that region and it shows that we have to consider the contribution of the 
zero mode. The zero mode means to consider the mode where the momentum, k, is 
zero. Therefore we have 



We can write the field for zero mode in terms of creation and annihilation operators 
as 




(3.4) 



X (t,t') = u a + u* al , 



(3.5) 



where 



[a , a ) = [oj, a$] = ; [a , aj] = 1 I [ a o, %^o\ = ■ 



(3.6) 
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Now we can calculate the commutator for the zero mode which is given by 



<0|[X (t),X (t')]|0> = |aVr: . (3.7) 



In the limit where /i — > we get 



(0\[X (t),X (t')]\0) = ^a'7rr* c At . (3.1 



Including the zero mode commutator the Green's function is given by 

G(r, t; r', t>) = -i ( (0| [X(t, r),X(t', r')] |0) + (Q\[X o (t),X o (f)]\0)) . (3.9) 



Thus we have 



/ : G (r,t;r',t') = —r:(r: + At) 



r 

/// : G(r,t;r',t') = 



// : G(r,t;r',t') = |p(r* Z + r^-At 2 ) 



IF : G(r,t;r',t') = ^(r* 2 + r*' 2 -At 2 y, 

V : G(r,t;r',t') = ^^r* c (r:-At) . (3.10) 

We can see that the Green's function satisfies the equation of motion in both vari- 
ables. We can also see that the retarded, t — t' > 0, and advanced, t — t' < 0, 
Green's functions appear with the opposite signs. Another observation is that the 
Green's functions in regions II and IV don't depend on the cut-off, r c , due to the 
contribution of the zero mode. We can also see that the Green's functions at regions 
I and V are zero when we set r c — > oo or r* — > 0. 

3.1.2 Hadamard Function 

The correlation function between two points in the space can also be explained using 
the Hadamard function of the field. For the mode expansion of the fluctuations of 
the string in AdS3 space, eq (2.30), the Hadamard function is given by 



(0\{X(t,r),X(t',r')}\0) = 2 dk 







Re(u k (t,r)ut(t',r')) (3.11) 



o 



+ / dkiu k {t,r)ul,{t',r'){0\al,a k \0) + u* k {t,r)u k it' ,r'){0\ala k ,\0) 
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As it's clear from its definition, for a free theory, the Hadamard function obeys the 
wave equation in both of its arguments (section 2.7 of [21]). According to (3.11) the 
Hadamard function depends on the state of the vacuum. AdS3 is not thermal and 
we have 

a*|0)=0. (3.12) 

So we get 



OG 



(0\{X(t,r),X(t',r')}\0) = 2 dk Re(u k (t,r)ut(f,r')) . (3.13) 

Jo 

Using the solution to the wave equation in AdS3 space, eq (2.28), one can easily 
see that if we consider massless modes where \x = 0, these integrals are divergent 
in the limit k — > 0. Therefore we should regularize the integrals. The logarithmic 
divergence in the integrals appears when k — > 0. Therefore we change the lower 
limit of the integral to fi ^ instead of zero. After the integration we can take the 
limit fi — > 0. The final result for the Hadamard function is 

(0\{X(t,r),X{t',r')}\0) = ^ ( - 8rf (In// + 7) + 12r* 2 - 4r*(r* + r*') - 4rV 

+ (At 2 - (r* 2 + r*' 2 )) \n\(At 2 - (r* - r*'f)\ 

- {At 2 - (r* 2 + r*' 2 )) ln|(-2r* + r* + r*') 2 - At 2 \ 

- 4r*(r* + At) ln| - 2r* c + r* + r*' - At\ 

- 4r*(r* - At) ln| - 2r* + r* + r*' + At\j . (3.14) 

We can simply see that the divergent term (ln/i) and the term coming from the 
regularization (7) are not coordinate dependent and can be thrown away. So we will 
be left with a finite result. 

In addition to the above result we also have to consider the contribution of the 
zero mode (3.4). The Hadamard function of the zero mode is 

(0\{X o (t),X (t')}\0) = ^< , (3.15) 



which in the limit \i — > reduces to 



—In' 1 

(0\{X (t),X (t')}\0) = -^rrr:- . (3.16) 



This should be added to our previous result for the anti-commutator. But you see 
the contribution of the zero mode is divergent but coordinate-independent. So it 
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can be thrown away. Therefore our final result for the regularized AdS3 Hadamard 
function which we denote as g(t,r;t' ,r') is 



g(t, r; if, r') = ^ ( Urf - 4r* (r* + r*') - Ar*r*' 



+ (At 2 - (r* 2 + r*' 2 )) ln|(At 2 - (r* - r*') 2 )\ 

- (At 2 - (r* 2 + r*' 2 )) ln|(-2r* + r* + r*') 2 - At 2 | 

- 4r*(r* + At) ln| - 2r* + r* + r*' - At\ 

- 4r*(r* -At) ln| - 2r* + r* + r*' + At\) . (3.17) 



We can see that it satisfies the equation of motion for the modes on AdS3 background 
as expected. 



3.2 BTZ Two-Point Functions 

In subsection (2.2) we explicitly wrote the solution to the wave equation on BTZ 
background. So we are now able to calculate the two-point function of the fields. As 
has been mentioned in the beginning of the section the commutator of the fields is 
state-independent and we can calculate it without being worried about the thermal 
vacuum of the BTZ space. But this is not the case for the Hadamard function and 
the effect of Hawking radiation should be considered. 



3.2.1 Green's Function 

In the previous section we already discussed the details of how to calculate the 
Green's function. In this section we skip the details and concentrate on the final 
results. We use the BTZ wave equation solution, eq (2.38), to write the commutator. 
One can see that the integrals are very complicated. So in order to be able to obtain 
the Green's function explicitly we consider the limit where r c — > 00. Therefore the 
commutator gets very simplified and reduces to 

{0 \[X( M r'Xm} = I dk - r , r(r / +tf / 4) (3.18) 

l k 2 l A coskr* coskr*' + rr' sinkr* sinkr*' 

+ kl 2 r coskr* sinkr*' + kl 2 r coskr*' sinfcr* 

The integrals are finite as it should be. Therefore the BTZ Green's function in the 
limit r* — > 0, in different regions of figure (2), gets the form 

I : G(r,r',At) = , 
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II 
III 

IV 

V 



G(r,r',At) 
G(r,r',At) 

G(r,r',At) 
G(r,r',At) 



2na'l 



H 



1 - cosh( — — ) sech( — — ) sech(— — 



r- 



-2-koII 



H2 V 



H 



1 — cosh( 



r H At, 



sech( 



TjjT 



sech( 







(3.19) 



We can see that the Green's function satisfies the equation of motion. Similar to 
AdS3 Green's function in r* — > limit, the Green's function is zero in regions I and 
V. It is also zero in region III which is outside the light cone. The retarded and 
advanced Green's functions also appear with opposite signs. 



3.2.2 Hadamard Function 

To calculate the Hadamard function of the transverse modes in BTZ background we 
use the same equation as in the AdS 3 space, (3.11), where the BTZ wave equation 
is given in (2.39). Note that in BTZ case the vacuum is thermal and we have 

< a l«*> = ' (3 - 20) 

where /3 — ^- and T H is the Hawking temperature of the black hole, T H = 
Therefore 

4r// 2 r°° 1 Bk 

(0\{X(t,r),X(t',r')}\0) = — jf dk ^ - fc2/4) cotlA cosfcAt (3.21) 

{kl 2 coskr* + r smkr*){kl 2 coskr*' + r' sinkr*') . 

In order to simplify the calculations we would like to obtain the Hadamard 
function in the limit where r — > oo. Therefore the integral can be expanded in 
terms of r* as 

Aa'l 2 f°° kB 
(0\{X(t,r),X(t',r')}\0) = -^- / dk coskAt coth^ (3.22) 

r H Jo 2 

^k(r 2 H + k 2 ty 3 r*> 3 - ^(4r% + hr 2 H l 4 k 2 + Z 8 fc 4 )(r* V' 5 + r* V' 3 ) + 0(r* 10 ) . 

Let's consider the first term in the expansion. The integral over k can be simplified 
to the form 

9(t, r; If, r') = -% r* V' 3 jf dk k(r 2 H + k 2 l*) , (3.23) 
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where we have used the following relation 



coth 



(3k 

T 



1 - e-? k l-e? k ' 



1 1 



(3.24) 



The result of the integral is 



g(t,r;t',r') 




sinh 4 ^ 



1 



(3.25) 



where we have used the fact that 8 = ?p- = It can be observed that the 

Hadamard function is periodic in imaginary time, which it should be since the 
vacuum is thermal. 

4 Two-Point Functions in Vaidya Background 

In the previous section we mentioned that the Hadamard function of quantum fields 
which is the expectation value of the anti-commutator depends on the choice of the 
vacuum. Therefore in a black hole spacetime it encodes the information about the 
Hawking radiation. In a system with a global time-like killing vector the vacuum 
is the zero eigenstate of all annihilation operators. The creation and annihilation 
operators are associated with positive and negative frequency modes with respect 
to the time coordinate of the system. In a Vaidya spacetime one can not define one 
global time-like killing vector for the whole space and therefore the choice of positive 
frequency modes and the definition of the vacuum is ambiguous. 

To overcome this problem one has to calculate the Bogoliubov coefficients which 
relate the complex basis of the fields which are solutions to the equation of motion 
in the vacuum spacetime to the basis in the black hole region. In general it is 
hard to do this calculation. Therefore people have done it using geometric optics 
approximation in the late time limit where the spacetime is stationary. This is the 
whole idea of the late time calculation of Hawking radiation. 

In the set-up of this paper we need to address the on-set of Hawking radiation 
which means to calculate the Hawking radiation in a finite time when the black hole 
forms. In the AdS/CFT language this calculation captures the h effects. One of 
the elegant ways of calculating the on-set of Hawking radiation is the work done by 
Callan, Giddings, Harvey and Strominger (CGHS) [18] in which they use the trace 
anomaly method. The set-up of their problem is the formation of a 2-dimensional 
black hole in linear dilaton background. This method is not applicable to the BTZ 
black hole formation due to the fact that the 2-dimensional worldsheet problem does 
not preserve diffeomorphism invariance 1 . 

lr This can be seen from the r 2 dependence in front of the third component of the 3-dimcnsional 
Vaidya metric which explains the excitations of the string. One can also write this system as a 
massive scalar field on the 2-dimensional Vaidya background. 
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In this section we introduce a fairly simple way to calculate the Hadamard func- 
tion in the Vaidya background which as far as we are aware, has not been used in the 
literature. The idea is to calculate the Hadamard function of the scalar field modes 
in the vacuum part of the space and then propagate it to the black hole region using 
the Green's function in the black hole space. The problem is an initial value problem 
where the value of the vacuum Hadamard function on the shock wave acts as the 
initial condition. We expect that this method would give us the information about 
the on-set of the Hawking radiation. We address this method as the propagation 
method in the rest of the paper. We will first explain the method more explicitly 
while applying it to the known example of CGHS [18]. Then we will use it in the 
BTZ black hole formation calculation. 



4.1 Linear Dilaton Black Hole Formation 

In this section we would like to apply the propagation method to the 2-dimensional 
set-up of formation of the black hole in linear-dilaton background from the collapse 
of a null shell [18]. The linear dilaton vacuum solution is given by 

~ 2 p = -\ 2 x + x~ , (4.1) 



e" 2 ^ = e 



where 

9+- = ^ , (4.2) 

is the metric in the conformal gauge and (x + ,x~) are the null coordinates. Using 
the following change of coordinates the metric reduces to the flat metric: 

x + = \e Xy+ ■ x- = -\e- Xy ~ . (4.3) 

A A 

The range of the coordinates x + and x~ are (0, oo) and (— oo, 0), respectively. There- 
fore y~ and y + both change over the range (— oo, oo). 

After the collapse of the null shell (shock wave) at the black hole 

solution is given by 

e -2<f> = e -2p = ax + _ X 2 x+x -> ^ ( 44) 

where a is the amplitude of the shock wave and 

x = x + — , (4.5) 

is the shift in the null coordinate x~ of the vacuum, due to formation of the horizon. 
Therefore the coordinate change to flat coordinates for the black hole is given by 

x + = \e Xa+ ; x~' = -\e- Xa ~ . (4.6) 

A A 

The relation between the coordinates in the vacuum and black hole regions is 

y + = o+ , 
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y- = -i ln(e-^- + . (4.7) 

As has been observed in [18] the difference between y~ and o~ coordinates is respon- 
sible for the nonzero and nontrivial energy-momentum tensor on which gives the 
Hawking radiation. 

The method used in [18] is based on trace anomaly argument 2 which gives a novel 
way of calculating a time-dependent Hawking radiation for the process of black hole 
formation. We would like to redo the calculation using the idea of propagating the 
two-point function (Hadamard function of the scalar field modes) from the vacuum 
to the black hole region. To do this we consider the value of the two-point function 
of the scalar field in the linear dilaton background, as the initial condition and 
propagate it to the black hole side using the black hole retarded Green's function. 

As the first step we need to obtain the Hadamard function and commutator of 
the scalar field in the linear dilaton and black hole spacetimes. We first need to 
solve the equation of motion for the scalar field. The wave equation of motion for 
the conformal matter is 

d + d^f(x + ,x~) = . (4.8) 
Therefore the solution can be decomposed into left moving and right moving parts 

f(x + ,x-) = f + (x + ) + f„(x-) . (4.9) 

Since in this framework the right movers transform the Hawking radiation we can 
ignore the left moving modes and write / in terms of creation/annihilation operators 
as 



f = duj [a u u u + alu*J , 

Jo 

poo 

= / du[b u v u + btvZ\, (4.10) 

J 

where and v u are the wave functions in the vacuum and black hole region, re- 
spectively and the creation/annihilation operators satisfy 

[a u ,a u >]=0; [a*, a*,] = ; [a m a*,] = 8(u - u') . (4.11) 

After imposing the boundary conditions at Z~[ and Zj£ the wave functions are given 

by 3 

v = 1 p-wiT 

1 e- iu}a ~ , (4.12) 



2 We will explain this method in more details in the appendix A. 

3 Note that for the mode expansion we have followed the notation of [22]. 
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where they satisfy the normalization conditions 

(u w ,u u >) = {v^Vu') = 2tt5(u - u') , 

= (uwX0 = O. (4.13) 

Using this mode expansion one can write the Green's function in linear dilaton 
background as 

Gv{yt,V2) = — « ( [/- (2/r ) > /- (2/2 )] > > 

= 2^ Im(u w (yf)u*(%)) , 

= si S n (2/i -2/2 ) ■ ( 4 - 14 ) 
Similarly for the black hole background we have 

_ 7T _ 

GbhK > CT 2 ) = - g si g n ( a i - ff 2 ) • ( 4 - 15 ) 
The anticommuator of the modes in the linear dilaton background is given by 

({/-(yD,/-(y a ")}> = a^ReWyrKCyi")). 

= -7 - In// - ln|yf - y^\ ■ (4.16) 

Note that the first two terms which come from the regularization are not y~ de- 
pendent and can be ignored. Therefore the regularized anticommutator is given 
by 

gvacuumiV 1^2) = ~MVl ~ Vll ■ ( 4 - 17 ) 

Now we have the necessary tools to calculate the Hadamard function after the 
formation of black hole using the propagation method. We would like to propagate 
each point of the vacuum Hadamard function to the black hole space. In general 
the propagation formula is 

(f)(u,v) = J du'dv' G(u,v;u',v')J(u',v') , (4.18) 

where (u, v) are null coordinates, G is the retarded Green's function of the space 
and J(u,v) is the source. Let's assume that the shock wave is nonzero at v = 0. In 
our set-up there is no source and for v > we have □</> = 0. But satisfies the 
initial condition on v = hypersurface (shock wave) which is (f)(u,v = 0) = <po{u). 
Since is also zero for v < we can assume that it satisfies the following equation 

□0 = d u Mu)5{v) . (4.19) 

Comparing this equation with the general Poisson's equation D(f>(u, v) = J(u, v) we 
can consider the right hand side of (4.19) to act as the source. 
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In the set-up of 2- dimensional linear dilaton black hole formation is in fact the 
black hole Hadamard function g and 0o is the vacuum Hadamard function calculated 
on the shock wave (x + = Xq). We can propagate each point from the shock wave 
using the retarded Green's function, separately. The equation of motion for the 
retarded Green's function in black hole space is 

d a+ d a -G(a + , a'; a + ' , a') = 5(a + - a + ) 8{a~ - a~) , (4.20) 

where 

G(a + , a'; a + \ a'') = --G(a~, a'') 0(a+ - a + ') Q{a' - a') . (4.21) 

7T 

The propagation formula (4.18) for the 2-dimensional set-up of this section reduces 
to 



fiV-bh = J dat'dcr-L ' G(crf, a 1 ; , a 1 ') d <T -g {y 1 '; y 2 ') 5{a^' - a%) , 



(4.22) 



where g w -bh means the Hadamard function with one point in the vacuum region and 
the other point in the black hole side. Due to the presence of the S(v) function in 
the definition of the source the integral over da + simplifies. We can also see that 
because of the simple form of the Green's function in black hole space (4.15) the 
integrand is a total derivative. Therefore the Hadamard function after the formation 
of the black hole is given by the vacuum Hadamard function calculated at the end 
points of the integral. Note that due to the presence of the function the range of 
the integration variables a/ is (— oo,o"j). 

In order to compare our result with the result in CGHS paper we need to calculate 
the energy-momentum tensor T which is defined as 

<0|T_(O|0> = <O|i0_/(O&-/(O|O> • (4-23) 

On the other hand one can see that T is equal to the regularized ^d a -d a -g(a^ , a 2 ) 

where a± — > a 2 . The derivative of the Hadamard function is 

d c- d o-9{Pii°2) = d a -d a -g(y^(a^),y 2 (a 2 )) , (4.24) 



<9 -<9 -In 



-1 

T 



ln(^ + e^)-ln(Y + e- A ^ 
A A 



-A 2 



-2 

ln(f + e- Xa i) -ln(f + e" 



;i + fe^r)(l + | e 



In order to obtain the energy-momentum tensor we have to set a l = a 2 . But it can 
be easily seen that the Hadamard function is divergent in this limit. To sort out this 
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problem we use point splitting method to eliminate the divergency. Let's assume 
a 2 = °\ + ^ an d then send 5 to zero. Therefore the above result is regularized 
by subtracting the term i. This corresponds to the contribution of the black hole 
vacuum [22] . Therefore we will be able to get the normal ordered energy momentum 
tensor (in leading order in 5) which is 



<T__) 
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1 



(4.25) 



This is exactly the result obtained in [18 
a 

Hawking radiation. 



Near the horizon which corresponds to 
— > oo the stress tensor reaches the constant value of ^ which is the late time 

4o 



4.2 BTZ Black Hole Formation 

As it was mentioned before we can not use the trace anomaly method to calculate 
Hawking radiation in our set-up which is the formation of the BTZ black hole from 
the collapse of a null shell. One could look at the fluctuations of the string on the 
other directions of the boundary for example one of the directions on the torus in 
the background AdS3 x T 4 x S 3 . This calculation is similar to the calculation for a 
conformal scalar on an AdS2 Vaidya background which has been done in appendix 
A. We will see the instantaneous thermalization of the system after the formation 
of the black hole which is expected from the discussions of [20] . 

In this section we would like to repeat the calculation done in the previous section 
to get the Hadamard function of the BTZ black hole on the boundary. In order to 
simplify the problem we consider the limit where r c — > oo. So we can compare the 
result with (3.25). This corresponds to an infinitely massive quark on the boundary. 
In this limit the scalar field X which describes the fluctuation of the string (quark 
on the boundary) behaves like r* 3 close to the boundary. This corresponds to the 
normalizable mode. 

Similar to the previous calculations we start from the AdS3 Hadamard function 
in AdSs background where one of the points is on v[ = 0, the shock wave, and 
the other point has the general value v 2 . The analogy is to propagate the v[ — 
point to BTZ boundary. Then we set v' 2 = and propagate the second point to 
the boundary. We will see that this sorts out the problem of divergency in AdS3 
Hadamard function when v[ = v' 2 = 0. 

Therefore the BTZ Green's function (G) and AdS3 Hadamard function (g) in 
the limit r c — > oo get simplified to 

G(t,r*;t',r*') = G(t, r*; t', r*') 0(u' - v)G{u - u')G(v - v') , (4.26) 
G(t,r*;t',r*') = (l - cosh(^f (t - t')) sech(^f r*) sech(^f r*')) , 



a' 



g(w 1 ,v 1 ;w 2 ,v 2 ) = 



- (wi - v 1 )(w 2 - v 2 )+ 
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(coxvt + u 2 v 2 - ^(wi + ^i)(w 2 + U2)) In 



(wi - w 2 )(v 1 - f 2 ) 



(w 2 - Vi)(Wi - f 2 ) 



where w = t — r* is in fact the coordinate u in AC1S3 space. The relation between w 
in ACIS3 and u in BTZ is 

9/2 „ 

w = _L tanh(-| M ) . (4.27) 
r# 2r 

In order to simplify the calculations we will work with the modified two-point func- 
tions defined as 

G'(t,r*;t',r*') = rr'G(t,r*;t',r*') , 

g'(t,r*;t',r*') = rr'g(t,r*;t',r*') . (4.28) 

Note that the modified BTZ Green's function satisfies the following equation of 
motion 

d v d u G' + ^ — ^r^ G' = 5{u - u')5(v - v') . (4.29) 

Therefore a general solution to the equation of motion, <j)(u,v), with the source 
J(u, v) satisfies the equation 

and can be written as 

4>(u,v) = J du'dv' G(u,v;u ,v')J(u' ,v') , (4-31) 

where G(u,v;u' ,v') is the Green's function satisying the equation (4.29). Similar 
to the discussion in subsection 4.1, one can assume that the source gets the form 
J(u,v) = d u <fto(u,v) S(v). Therefore <fi can be written in the form 

(f)(u,v) = I du'G{u,v,u',v' = 0) d u ,<f) {u') , (4.32) 



In the BTZ formation calculation <p is in fact AdS3 Hadamard function g' where 
v[ = and G is BTZ Green's function, G . Therefore we have 

POO 

#AdS 3 ;BTZ = / dtti 6(«i-Vi)6(«i-«i) 

Jo 

G'{ui, vi, u[, v[ = 0) d u > i g'(u' 1 ,v' 1 = 0; u 2 , v' 2 ) . (4.33) 

<?AdS3;BTZ is the two-point function where (ui,vx) is in BTZ space and (u 2 ,v 2 ) is in 
AdS3 space. The functions in the integral change the limits of the integration 
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to (v,u) which is (t + r*,t — r*). Therefore we can simplify the integral further by 
changing the integration variable as 

u[ = t x - x x r* . (4.34) 

Therefore we get 

#AdS 3 ;BTZ = -r\ J dxi G' (u x , v x ] u[, v[ = 0) d u ' i g\u' 1 ,v[ = 0;u 2 ,v' 2 )\ u > i=tl _ Xir * . 

(4.35) 

Since we are interested in what happens close to the boundary and we have already 
sent the boundary to infinity, we expand the integrands around r* — > 0. Therefore 
the integral is more simplified and the final result after integration is 

#AdS 3 ;BTZ = J^P? ^r H coth( r ^)d Ul g(u 1 (ui)]U2,v 2 ) 

+ l 2 d? 1 g(u) 1 (u 1 );u 2 ,v 2 )J\ ul=tl . (4.36) 

Note that one can take into account the fact that the derivatives of the g with 
respect to BTZ coordinate u should be changed to the derivatives with respect to 
AdS3 coordinate, u>. We can calculate gAds 3; BTZ explicitly. Its value at v' 2 = will 
be the initial data for the propagation of the second point to the boundary 

2a'r| r 2 uj\ 

SAdS 3 ;BTZ = Pi ( _ 2/2( _ 1 + cosh( ^)) + ^ S inh(^))2 ' 

where we have used (3.17). Now we can use the same integral (4.35) and propagate 
the second point to the boundary. The final result after propagation is 

/ _ a r H *2 *2 1 ( A oo\ 

9BTZ ~ 6I« ri ^ sinh 4 (^l) • ( 38) 

Therefore the BTZ Hadamard function is 

to = m r ^ s^hi^pij ' (439) 

One can see that this result is exactly the same as the result we got before for 
BTZ Hadamard function (3.25). As we can see here the BTZ Hadamard function is 
periodic in imaginary time with the period T# . Therefore it represents a thermal 
state at temperature Tjj. This implies instantaneous thermalization on the boundary 
of the 3-dimensional Vaidya spacetime for non-local operators like the Hadamard 
function. 

One could also assume a detector at fixed r which will immediately after the 
shock wave will feel the thermal bath. 

This result is also similar to the result obtained in [20] for the correlation function 
of the boundary operator for AdS 2 spacetime, using AdS/CFT duality. In the small 
separation limit where At <C 1 it behaves like ^ an d for large separation At ^> 1 
it decays exponentially e _4Ai . 
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Appendix 

A AdS2 Black Hole Formation 

In this appendix we show how the calculation for trace anomaly on AdS2 black hole 
formation predicts the rapid thermalization. This can be seen as the fluctuations of 
the string in the other directions of the space AdS,3 x M where M can be T A x S 3 . 
For example the transverse direction can be chosen to be along one of the torus co- 
ordinates. This problem is similar to solving a scalar field equation on AdS2 Vaidya 
background. The system satisfies all the necessary conditions to use trace anomaly 
method; there is a conformal scalar and the set-up preserves diffeomorphism invari- 
ance. We will do this calculation using both trace anomaly method and propagation 
of the Hadamard function of the scalar field from the vacuum to the black hole space. 
The energy-momentum tensor of the 2-dim theory is 

= (^-g^g p °G IJ d p X I d a X J + Gjjd^X'dvX^ , (A.l) 

where we have used the action (2.2). Therefore we explicitly get 

1 r 2 

T tt = ^^(dtXdtX - gld r Xd r X) , 

T rr = -^^{d r Xd r X - g 2 rr d t Xd t X) . (A.2) 

One can also write the energy-momentum tensor in null coordinates which would be 

T =0 

■ L uv u ) 

1 r 2 

Tvv = 2^P {9vXdvXh 



26 



T uu = ^^(duXd u X). (A.3) 

Using these coordinates one can see that the shock wave in the Penrose diagram in 
figure (l) 4 travelling along u direction with the amplitude a corresponds to 

T vv = a5(v) . (A.4) 

Using the argument of trace anomally the expectation value of the energy-momentum 
tensor is given by 

CO = £ ■ (A.5) 

where c is the central charge which is one for the conformal matter (scalar field) X 
and R is the ricci scalar which is the constant ^ for the AdS2 and asymptotically 
AdS2 black hole spacetime. We work in the conformal gauge where we can write 
the metric as ^ 

9vu = p ] 9w = 9uu = . (A-6) 
Therefore using (A.5) the conservation of the enrgy-momentum tensor implies 

(T vu ) = -T7—d v d u p , 

127T 

(T m ) = -—(d v pd v p-d 2 p + t v (v)), (A.7) 

1Z7T 

(T uu ) = -77 z-{d u pd u p- dlp + t u (u)) . 

1Z7T 

where t v and t u are integration constants which should be fixed by boundary con- 
ditions. Note that for the empty AdS2 space, where there is no nonzero energy flux 
or radiation, both T uu and T vv are zero while T uv is nonzero and produces the AdS2 
curvature. 



A. 0.1 Scalar Field on AdS2 Vaidya Background 

The first step to calculate the energy-momentum tensor is to find the right coordi- 
nates for the Vaidya background which keeps the metric continuous across the shock 
wave (metric itself, not necessarily its derivatives). We can use the calculation in 
section 2 which led to (2.10). Assuming the metric has the general form 

ds 2 = -e 2p dudv , (A.8) 

4 Note that the Penrose diagram for the AdS2 black hole formation is similar to the BTZ black 
hole formation except that there is no singularity in the AdS2 case and the black hole is explained 
by the coordinate which covers only a part of the space [20] . 
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in the two spaces of the problem we get 



ia * = -r *tf®(v - „» dvdu ' v<0 (A - 9) 

_4/2 

ds 2 = — 2 dvdu ; v > 

w cosh(^fu) - p| sinh(^|M) 

where u is the null coordinate in AdS2 black hole space. One can see that the metric 
is continuous across the shock wave v = 0, as it was expected. Therefore we get for 
v > 

(M 2 - dip = (d u p) 2 - dip = J , (A.10) 



and for v < 



(d v p) 2 -%p = 0, 



(M a -a2p = ^- (a- 11 ) 

At this stage we have to impose the boundary conditions to fix the integration 
constants t v and t u . We know that in the vacuum region the T uu and T vv are zero. 

r 2 

This implies that t v — and t u = — Therefore in the black hole region we get 



48tt/ 4 ' 

(T uu ) = . (A.12) 

This shows that the energy-momentum tensor component T vv is zero for v < but 
nonzero and constant after the shock wave. This happens due to the collapse of the 
null shell. Note that in contrast to linear dilaton black hole formation in [18], the 
energy-momentum observed in the black hole region changes from zero to a constant 
value instantly. Therefore one can conclude that the observer on the boundary 
instantly feels the thermal bath of the Hawking radiation due to the formation of 
the black hole. This is similar to the argument which has previously been given 
in [20]. 

We can also see that the curvature, R = 8e~ 2p 9 1 ,c? a ;p is constant and the same 
for all values of v and is given by 

In the next subsection we will compare the trace anomaly method result with 
the propagation of the Hadamard function. 



A. 0.2 AdS2 Vaidya Spacetime Two-Point Functions 

In order to use the propagation idea we first need to calculate the two-point functions 
in AdS 2 and AdS 2 black hole spaces. Since we just have the (t, r) coordinates of 
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AdS3 and BTZ spaces the equation of motion simplifies and the wave equation is 

1 2a' 

Uu(t,r*) = \ sinu;(r* — r*) e 



-ILut 



(A.13) 



The form of the wave equation is the same for both spaces and the only difference 
is in the definition of r*. So the Green's function (G) for both spaces is given by 



G(t,r*;t',r*') 



a ix 



sign(At + r* - r*') - sign(At - r* + r*') (A. 14) 



+ sign(At + r* + r*' - 2r*) + sign(At - r* - r*' + 2r*) 

where At = t — t'. One can see that the retarded Green's function which is 

G(v, u; v' , u) = G(v , u; v', u) Q(u — u) Q(v — v') , (A. 15) 

TlOi 

satisfies the equation of motion with the S function source 

d v d u G = 5(u - u) 8(v - v') . (A.16) 

The Green's function in different regions of the light cone diagram, figure 2, are 
given by 



I 

II 
III 

IV 

V 



G(t,t';r*y') = 0, 
G(t,t';r*,r*') = a'lr 

g(mW*') = o, 



(A.17) 



G(t,t';r*y')-- 
G(t,t';r*,r*') = . 



a tc 



The other two-point function is the Hadamard function of the scalar field in the 
AdS2 and AdS2 black hole space (similar to the Green's function the difference is 
only in the definition of r*) which is given by 



g{t x ,rl\t 2 ,r* 2 ) = a' 



- ln|ti - t 2 + r* - r* 2 \ - \n\ti — t 2 — r\ + r 2 \ 



(A.18) 



+ \n\tx -t 2 + rl+r* 2 - 2r*| + ln|^ - t 2 - r* x - r* 2 + 2r*\ 



which also satisfies the equation of motion for each variable. 
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Now we have necessary tools to calculate the Hadamard function of the scalar 
field in the AdS2 Vaydia space. We would like to propagate each point of the AdS2 
Hadamard function to the AdS2 black hole space. We will follow the same analogy 
as section 4. The initial condition in the propagation formulas (4.18) and (4.19) is 
given by AdS 2 Hadamard function calculated on the shock wave (v = 0). One can 
see that the AdS2 Hadamard function (A. 18) is divergent when one sets V\ — v 2 — 0. 
To avoid this divergency we propagate each point individually, which means we first 
set v[ = and propagate the first point and then we set v 2 — and propagate the 
second point. These words translated into formulas give 



9v-bh 



1 r°° 

— / du\ 9(«i - u' 1 )9(u , 1 - Vl + 2r*) 
ixa' J 

G(vi,ui, v[ = 0, ui) d ull g{v' x = 0, u[,v 2 , tu' 2 ) 



(A.19) 



Using the relation for the Green's function (A. 17) and considering the range of u\ 
imposed by the functions the above propagation relation reduces to 



flv-bh 



du[ d uli g(v[ = 0,u[,v' 2 ,u 2 



(A.20) 



where the relation between the null coordinate u in AdS2 black hole and the null 
coordinate u of AdS2 are given in (2.10). One can see that the integrand is the 
total derivative and the result is the AdS2 Hadamard function calculated in the end 
points of the integral. Note that the terms in g (A. 18) which don't depend on u[ 
cancel which means that the divergent term, ln|t>^ — v' 2 \, will be removed. Therefore 
we have 



9v-bh 



—a 



ln|a;i(«i) — ^1 — H^i^i — 2r*) — u' 2 \ 



\ia\u!i(ui) 



2r*\+ln\u 1 (v 1 -2r*) 



(A.21) 



Now we can propagate the second point. So we set v ' 2 = in g v _ bh and use the 
propagation formula (A.20) with the integration variable u 2 . One can see that only 
the terms which depend on u 2 will be left and we get 



^bh 



a 



lnjw^Mi) - u)2(u 2 )\ - ln\oJi(vi - 2r*) - u 2 (u 2 )\ 



ln|wi(iii) - ^2(^2 - 2r*)| + ln|wi(i;i - 2r*) - ^2(^2 - 2r*)| 



. (A.22) 



We are interested in the relation for the energy-momentum tensor. The energy- 
momentum tensor should be given by the regularized derivatives of the Hadamard 
function. In other words we have 



\Tu1u2) 



1 



Ana 



! du\ 9 1 regularized 



(A.23) 
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and the same for T VlV2 . The derivatives of the Hadamard function are 

Ft Pi - r ^ a ' 1 



4/ 4 sinh 2 ^^! -u 2 ) ' 

d Vl d V2 9 = tj^t] r- (A.24) 

4/ 4 smh - v 2 ) 

To regularize the above result we use the point-splitting method in which we assume 
U2 = U\ + e where e — > 0. The divergent part is which we can throw away. 
Therefore we are left with 

(T vv ) = (T uu ) = . (A.25) 

This is similar to what we got in the previous section using trace anomaly method 
(A.12). 
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